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Abstract 

We report on the results of a systematic study of X-ray flares from low-mass young 
stellar objects, using two deep exposure Chandra observations of the main region of 
the p Ophiuchi star-forming cloud. From 195 X-ray sources, including class I — III 
sources and some young brown dwarfs, we detected a total of 71 X-ray flares. Most of 
the flares have the typical profile of solar and stellar flares, fast rise and slow decay, 
while some bright flares show unusually long rise timescales. We derived the time- 
averaged temperature ((kT)), luminosity ((L x )), rise and decay timescales (r r and 
Td) of the flares, finding that (1) class I — II sources tend to have a high (kT), which 
sometimes exceeds 5 keV, (2) the distribution of (Lx) during flares is nearly the same 
for all classes from ~10 29 ' 5 to ~10 31,5 erg s -1 , although there is a marginal hint of a 
higher (Lx) distribution for class I than class II — III, and (3) positive and negative log- 
linear correlations are found between r r and Td, and (kT) and r r . In order to explain 
these relations, we used the framework of magnetic reconnection model with heat 
conduction and chromospheric evaporation to formulate the observational parameters 
(r r , Td, and (kT)) as a function of the pre-flare (coronal) electronic density (n c ), the 
half-length of the reconnected magnetic loop (L), and magnetic field strength (B). 
The observed correlations are well reproduced if loop lengths are nearly the same for 
all classes, regardless of the existence of an accretion disk. The estimated loop length 
is almost comparable to the typical stellar radius of these objects (10 10 -10 n cm), 
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which indicates that the observed flares are triggered by solar-type loops, rather than 
larger ones (~10 12 cm) connecting the star with its inner accretion disk. The higher 
(kT) observed for class I sources may be explained by a slightly higher magnetic field 
strength (^500 G) than for class II— III sources (200-300 G). 

Key words: ISM: clouds — ISM: individual (p Ophiuchi Cloud) — stars: flare 
— stars: pre-main-sequence — X-rays: stars 

1. Introduction 

Low-mass young stellar objects (YSOs) are classified into four evolutional stages based 
on the infrared (IR) to sub- millimeter spectral energy distributions (SEDs); the youngest and 
evolved protostars have class and I SEDs, while classical and weak-lined T Tauri stars (CTTSs 
and WTTSs) exhibit class II and III SEDs, respectively (Lada 1991; Andre et al. 1993; Andre, 
Montmerle 1994). In the 1980's, the Einstein satellite discovered X-rays from T Tauri stars 
(TTSs = CTTSs and WTTSs: Feigelson, DeCampli 1981; Feigelson, Kriss 1981; Montmerle 
et al. 1983). Koyama et al. (1994) first detected X-rays from the positions of class I sources in 
the p Ophiuchi dark cloud (WL 6 and WL 15 = Elias 29) with ASCA/GIS. Using ASCA/SIS, 
which has a better spatial resolution than GIS, Kamata et al. (1997) confirmed the X-ray 
emission from the class I sources. Independently, Casanova et al. (1995) obtained hints of 
X-ray detection from class I in the same star- forming area with ROSAT/PSPC. Then, succes- 
sive observations with ASCA further revealed that many class I protostars are X-ray emitters 
(Koyama et al. 1996; Ozawa et al. 1999; Tsuboi et al. 2000), thanks to a higher sensitivity 
for hard X-ray photons (>2 keV), which are less absorbed, because the extinction cross-section 
decreases as -Ex 2 ' 5 (Ex- X-ray energy). ROSAT observations also detected X-rays from class I 
(Grosso et al. 1997; Neuhauser, Preibisch 1997; Grosso 2001). They share the same character- 
istics of thermal emission with a plasma temperature of 0.5-5 keV and a strong variability with 
occasional rapid flares, consistent with the scenario of enhanced solar-type activity, attributable 
to magnetic dynamo processes. Furthermore, a recent observation with the Chandra satellite 
detected hard X-rays from class source candidates in the Orion Molecular Cloud 3 (Tsuboi 
et al. 2001), although its characteristics are still poorly understood due to the limited statistics. 
These pioneering discoveries demonstrate the unique capability of hard X-ray observations to 
probe YSOs or their close vicinity deeply embedded in dense cores. 

Based on the above results, the next step should be a systematic study of X-rays from 
YSOs and to approach physical conditions of these initial stages of stars, which gives extremely 
important information concerning star-formation theory. For example, our earlier results of a 
Chandra observation of the p Ophiuchi cloud suggested that class I protostars tend to show 
a higher plasma temperature (kT) than TTSs (Imanishi et al. 2001a, hereafter Paper I). 
Stelzer et al. (2000) also did a systematic study of TTS flares in Taurus-Auriga-Perseus with 
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ROSAT/PSPC, and found that the flare rate of CTTSs may be somewhat higher than that of 
WTTSs. To explain the quasi-periodic X-ray flare observed from the class I source YLW 15 in 
the p Ophiuchi dark cloud (Tsuboi et al. 2000), Montmerle et al. (2000) proposed a scenario of 
star-disk connection by one magnetic loop, where the central forming-star rotates faster than 
the inner edge of its accretion disk, which triggers periodic X-ray flaring. Shibata and Yokoyama 
(1999, 2002) further showed that kT is determined by the balance between reconnection heating 
and conduction cooling, and then derived the relation kT oc L 2 ^ 7 B 6 ^ 7 , where L and B are the 
half-length of the reconnected magnetic loop and strength of the magnetic field. Hence, larger 
L and/or B values make a higher temperature plasma. Combining these results, one plausible 
scenario is that younger sources have a much larger flare loop than evolved sources connecting 
the central source and the disk, and thus show frequent flares with a higher plasma temperature. 
Feigelson et al. (2002), on the other hand, reported that X-ray emission has no correlation with 
the presence or absence of the disk for the Orion Nebula Cluster X-ray sources. It is therefore 
still controversial whether there is any differences in the observed X-ray properties or emission 
mechanisms between the different classes. 

To address this issue, other important parameters would be the timescales of the flares, 
which were not considered in the above results. If we assume radiative cooling without suc- 
cessive heating (van den Oord et al. 1988), the decay timescale of flares (rd) is equal to the 
radiative loss timescale, defined as E^/R (E t ^ is the total thermal energy and R is emissivity 
of the plasma). Since E th is supplied by the magnetic energy (B 2 /8tt) and R depends on their 
plasma electronic density (Rocn 2 , n: electronic density), Td is tightly correlated with B and 
n. Also, a standard magnetic reconnection model (Petschek 1964) predicted that the rise 
timescale of flares (r r ) is proportional to the Alfven time (ta), defined as L/va, where va is the 
Alfven velocity (i>a oc B); hence, a larger L and/or a smaller B yield a longer r r . In fact, some 
stellar flares display unusually long rising phases (e.g., UX Ari: Giidel et al. 1999 and ROXs31: 
Imanishi et al. 2002). To examine any correlations between these observable parameters (kT, 
r r , and Td) would be fruitful for understanding of physical conditions of YSOs. However, little 
has been done for systematical approach so far. 

In the present work, we made the first systematic study of the flare activity of YSOs, 
using two deep exposure Chandra observations of the p Ophiuchi cloud (hereafter, p Oph) at 
a distance of 145 pc (de Zeeuw et al. 1999). Details of the observations are shown in section 
2. We then performed timing and spectral analyses for all detected sources, and examine any 
correlations between the derived parameters (section 3). From the observed quantities, we 
estimate the physical parameters and discuss the differences of the derived parameters between 
each class (section 4). 
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2. Chandra Observations of the p Ophiuchi Cloud 

The Chandra X-ray Observatory (Weisskopf et al. 2002) observed the central region of 
p Oph twice with a deep exposure of the ACIS-I array, consisting of four abutted X-ray CCDs. 
The first observation (here and after, obs.-BF) covered the south-east 17!4xl7!4 area, including 
cores B, C, E, and F, while the second observation (obs.-A) covered the north-west area centered 
on core A (Loren et al. 1990). Although some of the ACIS-S chips were simultaneously in 
operation, we do not use these data because large off-axis angles cause a degeneration of the 
sensitivity and position determination. From the Chandra X-ray Center (CXC) archive, we 
retrieved level-2 data, in which the data degradation caused by the increase of charge transfer 
inefficiency (CTI) in orbit was corrected. The X-ray events were selected with ASCA grades 
0, 2, 3, 4, and 6. The afterglow events were also removed. After processing, each observation 
yielded ^100 ks of live time (table 1). Earlier results of obs.-BF were found in Paper I. 

3. Analyses and Results 

3.1. X-ray Sources and NIR Counterparts 

Figure 1 gives an ACIS image of p Oph. The red and blue colors represent photons in 
the soft (0.5-2.0 keV) and hard (2.0-9.0 keV) X-ray bands, respectively. First, we discuss the 
source detection analysis using the wavdetect command (Freeman et al. 2002) in the CIAO 
package for 2048x2048 pixel images with the pixel size of ~0''5. We then consider the search 
for near-infrared (NIR) counterparts for all detected X-ray sources. The basic procedures of 
the analysis were the same in Paper I, except for the following points: 

(1) The significance criterion of wavdetect is relaxed from 10~ 7 (Paper I) to 10~ 6 . 

(2) We apply the wavdetect analysis for both the soft and hard X-ray bands as well as for the 
total (0.5-9.0 keV) X-ray band. 

(3) For reference of NIR sources, we use the Point Source Catalog in the 2MASS Second 
Incremental Data Release 1 , not that in Barsony et al. (1997), because the former has a better 
position accuracy (cr ~0'/l) than, and comparable sensitivity (Ks < 14.3 mag) to the latter 
(Barsony et al. 1997, a ~1"2 and K < 14.5 mag). 

We then detected 195 X-ray sources, nine of which were found in both the two obser- 
vations. Thirteen and nine sources were exclusively found in only the soft and hard bands, 
respectively (hereafter, the "soft-band sources" and "hard-band sources"). Table 2 (columns 
1-4) gives the name, background-subtracted ACIS-I counts, and the coordinates (after the offset 
correction) for each source. The soft-band and hard-band sources are indicated by the prefixes 
"S" and "H". The X-ray photons were extracted from a circle of l"2-23"3 radius, depending 
on the point spread function (PSF) radius for 1.49 keV photons, which is a function of the 

1 See ( http://www.ipac.caltech.edu/2mass/releases/second/doc/). 
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angular distance from the optical axis of the telescope (psfsize20010416.fits 2 in the CIAO 2.2.1 
package). Typically, in the source radius, ~95% of the X-ray photons are included. In order 
to obtain a high signal-to-noise ratio for some faint sources, we used a circle of a half radius of 
the PSF. Background counts were extracted from circles of 38 arcmin 2 and 57 arcmin 2 areas 
from the respective source-free regions in the ACIS-I fields of obs.-A and obs.-BF. 

In order to check for possible false sources with the wavdetect procedure (Feigelson et al. 
2002), we estimated the confidence level (CL) of the X-ray counts using Poisson statistics 
(Imanishi et al. 2001b). Then, 14 sources were found to have a CL value smaller than 99.9%. 
Although A-48 and A-H2 have significantly high CL, these may also be possible false sources 
because of the larger source size than the PSF radius (A-48) and of the severe contamination 
from A-2 (A-H2). We hence note that 16 sources were marginal detections, and label them 
with (m) in column 2 of table 2. 

About 60% of the X-ray sources had 2MASS NIR counterparts. Using these pairs, we 
determined and shifted the absolute positions of the Chandra sources so that the mean values 
of the Chandra-2MASS offset in the direction of right ascension and declination would become 
zero. The offset correction of obs.-BF (hence the position of the X-ray sources) was slightly 
different from those reported in Paper I because of the difference of the reference NIR catalogue. 
For the remaining sources, we further searched for counterparts using other NIR catalogues, 
and then identified an X-ray source, BF-90, to be a counterpart of GY 322 (Greene, Young 
1992). Finally, we concluded that 110 out of 195 X-ray sources have NIR counterparts. The 
offset between the X-ray and NIR sources is shown at column 2 in table 3. We also searched 
for radio (cm) and X-ray counterparts in published catalogues (columns 3-7 in table 3). 

Column 9 in table 3 gives the IR classification based on the spectral indices from the 
NIR to mid-IR band (ISOCAM survey at 6.7 and 14.3 /xm; Bontemps et al. 2001). Imanishi 
et al. (2001a) used the terminology class I c (class I candidate) for sources previously classified 
as class I (or flat spectrum source). In this paper, however, we regard all class I c s as class II, 
following Bontemps et al. (2001). Furthermore, we define some additional classes (class III C , 
BD, BD C , F, unclassified NIR sources, and unidentified sources). The definition of these classes 
is given in Bontemps et al. (2001), Imanishi et al. (2001b), and a footnote of table 3 of this 
paper. We then list 8 class Is, 58 class lis, 17 class Ills, 9 class III c s (class III candidates), 2 
BDs (brown dwarfs), 3 BD c s (BD candidates), 1 F (foreground star), and 12 unclassified NIR 
sources. Hereafter, we use the terminology "class III+III C " by combining the class III and class 
III C for brevity. 

2 This file was constructed by using the SAOsac raytrace code (Jerius ct al. 1995, see also 
http:/ /asc. harvard.edu/chart/). 
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3.2. Timing Analysis 

We made two X-ray light curves in the 0.5-9.0 keV region for all of the Chandra sources 
(background was not subtracted), with respective time bins of 2000 s and 4000 s. The source 
regions were basically the same as those used in an estimation of the source counts (subsection 
3.1). For the brightest two sources (A- 2 and BF-64), however, we used an annulus of 2."5-12."5 
and 2."5-7."5 radius in order to avoid photon pileup (Davis 2001). The light curves show many 
flare-like events. We defined a "flare time bin" by the following criterion: 

^>2, (1) 
AN ~ K 1 

where N, AN, and N are the X-ray counts, the statistical 1-a error for the relevant time bin, 
and average counts in a 20000 s interval including the relevant time bin, respectively. If the 
above criterion was satisfied in both of the two light curves (2000 s and 4000 s time bin), we 
defined the event to be a "flare". We then picked up 71 flares (table 4). Figure 2 shows all of 
the light curves of the detected flares. We note that the flare list in Paper I is slightly different 
from that in this paper, due solely to the different definition of a flare. This slight difference, 
however, has no significant effect on the following analyses and discussions. 

Most flares had the typical profile of those from low-mass main-sequence stars or YSOs; 
fast-rise and slow-decay, while some sources show unusual flares having slow rise timescale; 
e.g., BF-64 = YLW 16A (figure 2). We fit the rise and decay phases of the flare by a simple 
(exponential + constant) model using the QDP command in the LHEASOFT 5.0 3 , and derived 
the respective e-folding times and quiescent level (r r , Td, and Q; table 4) with their 90% errors, 
where the flare peak time (t p ) was fixed to be the maximum time bin, except for the second 
flare of BF-64. Evidently, this simple model can not reproduce the unusually giant flare from 
BF-64, even if we relax t p to be free. However, we do not intend to make a model of the flare 
light curve and only have interest in the typical timescales in this paper; hence, we do not 
discuss this discrepancy in further detail. 

3.3. Spectral Analysis 

X-ray spectra were made for all of the Chandra sources. For flaring sources, we made 
time-averaged spectra in the quiescent and flare phases separately. The flare phase was the time 
from {t p — r r ) to (t p + Td), and the quiescent phase was that outside of the flare phase. The 
background regions were the same as used in subsection 3.1. We then fit the spectra by a thin- 
thermal plasma model (MEKAL: Mewe et al. 1985) with the photoelectric absorption (WABS). 
The metal abundances were fixed to be 0.3 solar based on previous fitting results (Paper I), 
unless otherwise noted. If the temperature was not constrained, which is often the case for very 
faint sources, we fixed the temperature at two representative values of 1 keV (typical of TTSs) 
and 5 keV (protostars), then estimated the respective absorptions and luminosities. For the 

3 See ( http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/). 

6 



nine sources detected in both of the observations, we derived the parameters of each spectrum 
assuming the same absorption. This simple model was generally acceptable. BF-46 (R0Xs21) 
and BF-96 (R0Xs31), however, needed multi-temperature spectra with unusual abundances 
(Imanishi et al. 2002). A-23, A-24 and BF-10, on the other hand, required overabundances 
(subsubsection 4.8.1 in Paper I; Imanishi 2003). In table 2 (columns 5-9), we summarize the 
best-fit parameters; time-averaged temperature ((kT)), emission measure ((EM)), absorption 
column (iVn), flux, and luminosity ((Lx)) in 0.5-9.0 keV. 

3.4- Luminosity Function 

In order to examine the differences along the evolutional stages, we calculated the X-ray 
luminosity function for the detected X-ray sources of class I, II, and III+III C (figure 3), using 
the ASURV statistical software package (rev.1.2) 4 based on the maximum likelihood Kaplan- 
Meier estimator. The luminosity function of class I seems to be shifted toward higher values 
than those of class II and III+III C both in the quiescent and flare phases. 

To be more quantitative, we estimated the significance level using two nonparametric 
two-sample tests in ASURV: the Gehan's generalized Wilcoxon test (GW) and the logrank test 
(Feigelson, Nelson 1985). The results are given in table 5, where the blank (...) indicates 
that the significance level is less than 90%. From table 5, we can see that both the GW and 
logrank tests for the (Lx) difference between class I and the others in the flares show a marginal 
significance level of ~94%. We should further note that the effect of "undetected" faint flares 
of class I can not be ignored; if we assume an undetected class I has (Lx) smaller than 10 30 
erg s _1 , which is near to the detection threshold of a class I flare, and add this source to the 
GW and logrank test sample, then the relevant significance level is largely reduced to < 90%. 
We hence conclude that the higher (Lx) of class I than class II+III during flares is marginal. 

On the other hand, from table 5, we can see no significant difference between class II 
and III+III C source in the flare phase. This is contrary to the results of the Taurus-Auriga- 
Perseus samples (Stelzer et al. 2000), for which there was a significant difference between the 
distribution of flare (Lx) of class II (CTTS) and class III+III C (WTTS) sources. 

Also, we can see no significant difference of (Lx) in the quiescent phase among all classes. 
This is consistent with the previous estimation for this region with ROSAT derived by Grosso 
et al. (2000). They tried to estimate the significances rather strictly by considering the upper 
limit of undetected sources. Our samples give a more severe constraint because the detection 
threshold is largely reduced to give mean luminosities of io 29 - 5 -10 29,8 erg s -1 , which are about 
10-times lower than that of ROSAT. 

4 See (http://www.astro.psu.edu/statcodes). 
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3.5. Plasma Temperature and Flare Timescales 

In figures 4 and 5, histograms of (kT), r r , and t& are shown for each class separately, 
together with the mean values and standard deviations. We exclude the samples whose best- 
fit parameters and/or those errors were not determined, due to the limited statistics. For 
sources with multi-temperature plasma (BF-46 and BF-96), we regard (kT) of the soft and 
hard components as the quiescent and flare values, respectively, following the discussion that 
the former would be steady coronal emission, while the latter is the flare activity (Imanishi 
et al. 2002). 

Most of the (kT) values in the quiescent phase are in the range of 0.2-5 keV, while it 
becomes systematically higher (1-10 keV) in the flare phase, indicating plasma heating during 
the flare. Although the flare temperature of all classes is distributed around 3-4 keV, some 
flares of class I and II sources show a higher temperature than 5 keV, while all flares of class 
III+III C sources have (kT) less than 5 keV. The mean values of (kT) during flares for class I 
and II (4-5 keV) are therefore larger than that of class III+III C (2.7 keV). 

Like the case of Lx, we estimated the significance level for the (kT) difference among all 
of the classes with two nonparametric two-sample tests in ASURV. The results are also given 
in table 5. From table 5, we can see that both the GW and the logrank tests for the difference 
of (kT) during flares between class I+II and class III+III C sources show high significance level 
of 94% and 98%. One may be concerned that a systematically large absorption of class I 
and II sources makes artificial spectral hardening of these classes. We, however, argue that 
the observed tendency of higher (kT) during flares for younger sources is not changed by the 
extinction differences for the following reasons: 

(1) If the intrinsic distribution of the flare (kT) in class III+III C is the same as classes I and 
II, some of the class III+III C flares should show a higher (kT) than 5 keV, but there is no such 
flare in class III+III C . 

(2) We re-estimated the mean (kT) of class I+II and class III+III C during the flare for the 
limited sources of < 5xl0 22 cm -2 (most of class III+IIL, sources have lower A^h than this 
value). Although the mean Ah is nearly the same (2.5±0.2 and 2.2±0.3 x 10 22 cm -2 for class 
I+II and class III+III C ), the mean (kT) is still higher for class I+II (3.7+0.4 keV) than for 
class III+III C (2.4+0.3 keV). 

In the quiescent phase, on the other hand, ASURV shows a significance level >95% that 
the tendency of a higher (kT) of class I than the others (table 5). However, we suspect that this 
may be a bias effect because the quiescent temperatures are systematically low (< 5 keV), and 
hence would be significantly affected by the absorption compared with the flare temperatures. 

The rise and decay timescales of the flares (r r and ra) are distributed around 10 3 ' 5 s 
and 10 4 s, respectively. Although the mean value of r r indicates that younger sources have 
a shorter rise timescale, the difference of r r among these classes is not statistically significant 
(<90%, table 5). The larger mean values for class III+III C sources are primarily due to flares 
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with unusually long timescales of ^50 ks (A-2, A-63, and BF-96); such flares are not seen in 
class I and II sources (~30 ks at the maximum). Since we can not reject the possibility that 
these long timescales are simply because of the composition of two (or more) unresolved flares, 
we further estimate the mean r r and ra values of class III+III C sources without these flares to 
be 3.1±0.6 ks and 6.9±0.8 ks, which is comparable to those of the other classes. 

3.6. Correlation between the Flare Parameters 

Figure 6 shows the correlations between the derived parameters of flares; r r vs and 
(kT) vs r r . In these figures, we exclude those flares for which we could not determine the 
errors mainly due to limited statistics. Possible positive and negative log-linear correlations 
can be seen in r r vs and (kT) vs r r , respectively. We hence checked the significances of these 
correlations using Cox's proportional hazard model (Isobe et al. 1986) in ASURV. This gives 
significances of ~99% (r r vs T&) and ~94% ((kT) vs r r ). We further estimated the best-fit 
log-linear models with ASURV to be 

(Zl) = 10 2.5±0.4(Z£)0.4±0.1 (2) 

s s 

and 



(Zi) = 10 3.8±0.1 ( i^-0.4±0.3 (3) 

These models are shown by the solid lines in figure 6. 
4. Discussion 

From the detected 71 X-ray flares of the p Oph X-ray sources, we found that (1) the 
flares from class I and II sources tend to show a higher plasma temperature, which sometimes 
exceeds 5 keV, than class III+III c s; (2) the distribution of (Lx) is nearly the same among the 
classes, although there is the marginal hint of slightly higher (Lx) for class I sources; and (3) 
the plots of r r vs show positive and those of (kT) vs r r show negative log-linear correlations. 
Using these enormous samples and derived properties, we discuss the overall feature of YSO 
flares in this section. 

4-1. Flare Rate 

First, we estimate the flare rate (F) for each class. Following Stelzer et al. (2000), F is 
defined as 



F _ (f T + t d )N F Jaj+^y/Np 

NT obs NT ohs ' 1 ] 

where (f r , fd) and (cx Tr , a Td ) are the mean value and uncertainty of (r r , r<i) derived with ASURV 
(figure 5), N and N-p are the number of sources and detected flares, and T Q b s is the duration 
of the observations («100 ks). We then determine F to be 15±1%, 8.0±0.2%, and 13±1% 
for class I, class II, and class III+III C sources, respectively. All of the values are much higher 
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than that derived with ROSAT (F ~ 1%: Stelzer et al. 2000), which may be primarily due to 
the extended sensitivity in the hard X-ray band of the Chandra observations, because the flare 
activity (flux increase) is clearer in the harder X-ray band. 

The surprisingly large value of F for class III+III C sources is simply due to the effect 
of flares with unusually long timescales (A-2, A-63, and BF-96). In fact, if we regard those 
unusual flares as being affected by other unknown components (unresolved flares, for example), 
and exclude them, F of class III+III C becomes 5.0±0.1%, which is significantly smaller than 
those of class I — II sources. Such a tendency is consistent with the ROSAT results (F = 1.1±0.4 
% and 0.7±0.2% for class II and III). We hence suggest that younger sources tend to show 
frequent flare activity. 

4-2. Interpretation of the r r vs Td and (kT) vs r T correlations 

In order to explain the observed features of the flare parameters, we formulated r r , 
Td, and (kT) as a function of pre-flare (coronal) electronic density (n c ), half-length of the 
reconnected magnetic loop (L) and strength of magnetic field (B), based on the idea of a 
standard magnetic reconnection model (Petschek 1964) and the balance between reconnection 
heating and conductive cooling (Shibata, Yokoyama 2002). For simplicity, we assumed that 
is equal to the radiative loss timescale (r ra( j), although it would be possible that Td is slightly 
larger than r ra d (see appendix 3). We also assumed that the reconnection rate, Ma [ratio of the 
reconnection interval and Alfven timescales; equation (A4)], is 0.01, which is the mean value of 
Ma for the solar observations (0.001-0.1, table 2 in Isobe et al. 2002). Details of the formulae 
are given in appendix 1. 

4-2.1. r r vs Td 

From equation (A9), the reconnection model predicts a positive correlation between r r 
and Td as Td = At]I 2 . The slope of 1/2 shows good agreement with the observed value of 
0.4±0.1 [equation (2)], and hence supports our assumption that the decay phase is dominanted 
by radiative cooling. From equation (A9), we re-estimated log(A) with a fixed slope of 1/2 to 
be 2.06±0.04 (the dashed line in figure 6a). We also calculated log(A) separately for each class, 
and then obtain 2.1±0.1, 2.0±0.1, and 2.0±0.1 for class I, II, and III+III C , respectively. We 
can see no significant difference between the classes; hence, A is assumed to be the same for all 
sources. Since A is a function of n c [equation (A10)], we can determine n c to be 

iQlCUSiO.OQ/^A) [cm -3]_ (5) 

v o.or L J y ' 

Considering the possible range of Ma (0.001-0.1), this favors a higher pre-flare density for YSOs 
(10 9 -10 12 cm" 3 ) than that for the sun (~10 9 cm" 3 ). In fact, Kastner et al. (2002) derived the 
plasma density of a CTTS (TW Hydrae) using the ratio of Nerx and Oyn triplets, and found 
an extremely high density of ~10 13 cm -3 , even in their quiescent phase. Hence, a larger n c 
value for YSOs than that for the sun would be a common feature. 
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4-2.2. kT vs r r 

The predicted correlations between (kT) and r r are shown in equations (All) and (A12). 
Positive (r r oc (kT) 7 / 2 ) and negative (r r oc (kT)~ 7 / 6 ) correlations are expected for constant B 
and L values, respectively. The dash-dotted and dashed lines in figure 6b show constant B 
and L lines, in which we uniformly assume n c = 10 10,48 cm -3 [equation (5)]. Our observed 
negative correlation therefore predicts that flares with a higher temperature are mainly due to 
a larger magnetic field. However, the observed slope of — 0.4±0.3 is significantly flatter than 
—7/6 [equation (A12)], and hence the effect of L would not be negligible. We estimated the 
mean B and L for each class (table 6) using the mean values of (kT) and r r (figures 4b and 
5a), which are summarized in table 6. Younger sources tend to show larger B values; the mean 
values of B are ~500, 300, and 200 G for class I, II, and III+III C sources, respectively. This is 
consistent with the observed result of a higher plasma temperature for younger sources. Also, 
the estimated values of L (10 10 -10 n cm) suggest that L is nearly the same among these classes, 
comparable to the typical radius of YSOs. These results indicate that the flare loops for all 
classes are localized at the stellar surface, which is in sharp contrast with the idea of the flares 
for class I sources triggered by larger flare loops connecting the star and the disk (Montmerle 
et al. 2000). 

It is conceivable that the class III flares with unusually long timescales (A-2, A-63, and 
BF-96) mainly affect the systematically lower B for class III sources. Since it is possible that 
two or more flares make such unusually long timescales, we re-estimated B and L using the 
mean value of r r without these flares (subsection 3.5). The results are shown by the parentheses 
in table 6. Again, we confirmed the same results: class III+III C sources have a lower B than, 
and comparable L (but slightly lower) to class I and II sources. 

4-3. Comparison with the kT-EM Scaling Law 

Shibata and Yokoyama (2002) showed that L and B of flares can be estimated from the 
kT-EM plot (see appendix 2). This method is completely independent of the flare timescales (r r 
and Td), and hence can be used as a consistency check of our estimation derived in subsection 4.2. 
The (kT)-(EM) relation of the p Oph flares is shown in figure 7. Using ASURV, we determined 
the mean values of log((£M)) to be 53.77±0.20, 53.33±0.07, and 53.48±0.16 for class I, II, and 
III+III C sources, respectively. We then calculated the mean values of the magnetic field strength 
and loop length (-Bsy and Lsy) by equations (A15) and (A16) using the estimated value of n c 
= io 10 - 48±0 - 09 cm -3 (subsubsection 4.2.1). The results are given in table 6. The values of B$y 
(100-1000 G) are much higher than those derived in Shibata and Yokoyama (2002) and Paper 
I (50-150 G). This discrepancy is caused by two different assumptions concerning the n c value: 
10 9 cm -3 for Shibata and Yokoyama (2002) and Paper I (typical value of the solar corona), 
and 10 10,48 cm -3 for this paper. B^y and Lsy show good agreement with those derived in 
subsubsection 4.2.2 (B and L in table 6); higher _B SY values for younger sources and the same 
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order of L$y comparable to the typical stellar radius. Therefore, our estimation of n c , L, and 
B using the r r -Td and (kT)-r r relations is cross-checked. 

The same conclusion of subsubsection 4.2.2 and this section also indicates the propriety 
of our assumption of Ma = 0.01. If we assume a smaller (or larger) value of Ma, n c [equation 
(5)] becomes smaller (larger), which causes smaller (larger) estimations of -Bsy [equation (A15)], 
and makes a larger discrepancy between B and Bsy (table 6). 

4-4- Evolution of YSO Flares 

Combining all of the results discussed in the previous sections, we propose a simple view 
of the evolution of flare activity on low- mass objects as follows. In their earlier stage (class I), 
sources have a relatively strong magnetic field (?^500 G), and show frequent X-ray flares with 
a higher temperature (~5 keV). As stars evolve (class II and Ills), the magnetic field gradually 
decreases (200-300 G) and makes a moderate temperature plasma (2-4 keV) via X-ray flares. 
During these phases (class I — III) , the length of the flare loop does not change significantly 
(10 10 -10 n cm). As approaching to the main-sequence stage, the magnetic field and length of 
the flare loop become weak (50-150 G for the sun) and short (10 8 -10 9 cm), which causes a 
lower plasma temperature (0.1-1 keV) and shorter flare timescales (10-100 s), as can be seen 
in the sun. 

Generally, RS CVn systems show smaller flare timescales than YSOs (~100 s), while 
kT is comparable (~5 keV). Based on the above idea, this may be due to a larger B value. In 
fact, Donati et al. (1990) estimated B of RS CVns to be ^1000 G, which is as large as the 
maximum value for YSOs (figures 6b and 7). 

4-5. Comment on the Giant Flares 

Figure 7 indicates that the giant flares (EM ^10 55 cm -3 ) previously detected with 
ASCA such as V773 Tau (class III: Tsuboi et al. 1998) and ROXs31 (class III: Imanishi et al. 
2002), as well as BF-64 = YLW 16A in this paper (class I; EM ~10 55 cm -3 ), should have either 
a large L or B, if the sustained heating is negligible. Since the flares of ROXs31 and YLW 16A 
show a relatively large rise timescale (~10 4 s), they may have large L values (10 n -10 12 cm). 
Because the rise time of V773 Tau, on the other hand, is among the shortest (~10 3 s), the large 
Lx would be primarily due to the large B value. In fact, a detailed flare decay time analysis for 
V773 Tau predicted an extremely large B value of >1000 G (Favata et al. 2001). We suspect 
that the binary nature of V773 Tau (Welty 1995) causes such an exceptionally large B, which 
would be the same case as that for the main-sequence RS CVn systems. 

5. Summary 

We summarize the main results of the flare analysis of two Chandra observations of p 
Oph with fslOO ks exposure as follows: 
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1. From 195 X-ray sources detected in the main region of p Oph region, we found 71 X-ray 
flares. Most of the flares show the typical profile of the solar and stellar flares, while some 
bright flares have an unusually long rise timescale. 

2. Flares from class I and II sources tend to have a high plasma temperature, which sometimes 
exceeds 5 keV. 

3. There is a positive correlation between r r and Td, which is well explained by the standard 
magnetic reconnection model. Shorter r r and Td are due to the smaller L and/or larger B 
values. 

4. We found a negative correlation between kT and r r , which indicates the same order of the 
flare loop length regardless of their classes. Larger kT values for class I and II sources are 
due to larger B values. 

5. The expected loop length is comparable to the stellar size, indicating that the scenario of 
the star-disk arcade magnetic loop is unlikely. 
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Appendix 1. Magnetic Reconnection Model 

In this appendix, we show details for estimating the flare parameters (r r , Td, and kT), 
which are used in section 4. These formulations are based on the standard magnetic reconnec- 
tion model (Petschek 1964), as well as the observed results and analysis of solar flares. 

A. 1.1. Estimation of the Flare Parameters 

• Plasma temperature ((kT)) 
Shibata and Yokoyama (1999; 2002) showed from MHD simulations that the balance 
between the heating rate and conduction cooling determines the maximum temperature 
of reconnection heated plasma (T max ), and derived a relation for the temperature of the 
evaporated plasma (T) filling the reconnected magnetic loop: 

T ^l T ^ 6 , 5xm ^ r ^^ r{ _£_ r [K ], (A!) 

where n c , L, and B are the pre-flare (coronal) density, semi-length of the loop, and mag- 
netic field strength, respectively. As for the usual stellar flares, we could only determine 



13 



the time-averaged temperature ((kT)) due to the limited statistics. However, van den 
Oord et al. (1988) showed that the behavior of the flare temperature is exponential, hence 
we obtain the relation 

r 7^if^ ,Ba '' lff ' (A2) 

Using equations (Al) and (A2), (kT) is determined as 

( kT ) = 3 - 64 (T?/^)" 1/7 (T7iT^— ^ Wr ) 6/? t keV ]- ^ A3 ) 

10 y cm 6 10 11 cm 100 G 

Rise timescale (r r ) 

Based on a standard reconnection model, r r is equal to the interval of magnetic recon- 
nection. Petschek (1964) showed that the reconnection timescale is proportional to the 
Alfven time, ta = L/v\, where va is the Alfven velocity. Hence r r is 



_ r A _ ^Arrmn c L 
r ~M A ~ M A B ' 1 ' 

where m is the proton mass (=1.67xl0 -24 g). The correction factor, Ma, sometimes 
referred to as the reconnection rate, was estimated to be 0.01-0.1 (Petschek 1964), while 
observations for the solar flares show that Ma is in the range of 0.001-0.1 (table 2 in Isobe 
et al. 2002), regardless of their flare size. r r is therefore 

r, * 1.45 x 1 0^)-.,^)V 2( _A_)(_f_ ) - 1 W . (A5) 

• Decay timescale (ra) 

We assume that is nearly the same as the radiative cooling timescale (r ra( j), i.e, 

3nkT 

T ^ T ^^mTy (A6) 

where n and A(T) are the maximum plasma density and radiative loss function given by 
10 - 24.73 T i/4 for T > 20 MK (Mewe et al. 1985). Furthermore, we assume that magnetic 
pressure is comparable to the plasma pressure at the flare peak, 
B 2 

2nkT = — . (A7) 
8tt V ; 

Using equations (Al), (A6), and (A7), we obtain r d as 

A. 1.2. Predicted Correlations between the Flare Parameters 

• r r vs Td 

Using equations (A5) and (A8), we obtain a relation 



where 



-) = ^4(-) 1/2 , (A9) 
s s 
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Hence, a positive correlation is expected between r r and r d . 
• (kT) vs r r 
From equations (A3), (A5) and (5), we obtain 

^ 1 ° 2 ' 2 °'w''ei)" I (i#G)"'i§' 7/2 W. <A") 

* a 104 82 < w' I/3 <ll'"(lo^) 1/3 (l§)" /6 W- < A12 > 

These equations indicate that the positive and negative log-linear correlations of (kT) vs 
r r are expected if B and L are equal for all sources, respectively. 

Appendix 2. The kT-EM Scaling Law 

Using equations (Al), (A7), and assuming V = L 3 (V: plasma volume), Shibata and 
Yokoyama (2002) derived the relation between the flare maximum temperature (T) and emission 
measure [EM, equations (5)-(6) in Shibata and Yokoyama (2002)]. Similar to equation (A2), 
the time behavior of EM is also exponential (van den Oord et al. 1988), hence we assume 
EM = l.Q(EM). We thus replace equations (5) and (6) in Shibata and Yokoyama (2002) as 

w* 3 - 1 * i048 <i^> 3/2 O- 5 O 17/2 t-- 3 i' < Ai3 > 

(EM) = 7.0 x 10 5 ^-/^) 2 / 3 (-^) 5 /3(M)8/s [cm -3]. (A14) 
10 y cm a 10 11 cm kev 
We also derive the following two equations: 

B = 52(-^) 3 / 10 (-p^)-/ 5 (^) 17 / 10 [G], (A15) 
10 y cm d 10 &u cm ^ kev 

L - 7.8 x lO^--/^)" 2 /^--!^) 3 /^^)- 8 / 5 [cm]. (A16) 
10 y cm 6 10 ou cm 6 kev 

Appendix 3. Possible error for the derived parameters 

In the discussions of the main text (section 4), the assumption of [equation (A6)] 
includes relatively large uncertainty. Reale et al. (1997) proposed that r d becomes about a 
factor of < 10 larger than the radiative loss timescale (r ra d), if the sustained heating exists 
during the decay. The dependences of the relevant parameters on are: 

n c cx r d " 2 , (A17) 

B oc nl /4 oc r d " 1/2 , (A18) 

Locn~ 1/2 ocr]. (A19) 

Accordingly, n c has large dependence on r d ; larger n c value than ~10 10 ' 5 cm -3 may be conceiv- 
able. The dependence of L is also relatively large. However, this uncertainty makes L much 
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smaller, and hence larger flare loops are still less possible. 

Another uncertainty is in the kT-EM scaling law. Shibata and Yokoyama (2002) showed 
that the effect of the filling factor / (V = fL 3 ) may not be negligible. This gives possible errors 
for the values in columns 5-6 of table 6. The dependences of the relevant parameters on / are 



Hence, slightly smaller and larger values of .Bsy and L$y would be conceivable, although the 
effect is only a factor of < 5 (if / = 0.1). 
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Fig. 1. False-color ACIS image of p Oph. Red and blue colors represent photons in the soft (0.5-2.0 kcV) 
and hard (2.0-9.0 keV) bands, respectively. 
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Fig. 2. Light curves of all detected flares in 0.5-9.0 keV, binned in 2000 s interval. Names and classes for 
each flare are shown at the upper-left of the panels. The time axis starts at MJD = 51679.9944 (obs.-A) 
and 51647.7848 (obs.-BF). The solid line is the best-fit (exponential + constant) model. Since the light 
curve of BF-64 suffers photon pilcup during the second flare, we show the two flares separately (the former 
is extracted from the 7"5 radius and the latter is from 2"5-7"5). That of A-2 is also extracted from the 
2."5-12"5 radius circle. on 
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Fig.2 (Continued) 
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BF-78 (class II) 
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Fig.2 (Continued) 
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Fig. 3. Normalized X-ray luminosity functions of class I (dashed), class II (dash-dotted), and class III 
(solid) sources in the (a) quiescent and (b) flare phases. The mean value (Mean) and standard deviation 
(a) of logfix] m the un it of erg s _1 for each class are shown in the figures. The parentheses indicate errors 
of the mean values. 
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Fig. 4. Histograms of kT in the (a) quiescent and (b) flare phases for each class, with their mean values 
(Mean) and standard deviation (a) in the unit of keV for each class. The parentheses indicate errors of 
the mean values. 
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Fig. 5. Same as figure 4, but for the flare rise and decay timescales, with their mean values and standard 
deviation in the unit of ks for each class. 
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(a) T r vs T d 



(b) <kT> vs T r 




Fig. 6. Relations between (a) r r vs Td and (b) kT vs r r . The circles, triangles, squares, diamonds, and 
crosses represent flares from class I, class II, class III+III C , unclassified NIR sources, and unidentified 
sources, respectively. The solid lines represent the best-fit log-linear correlations derived with ASURV 
[equations (2) and (3)]. The dashed line in (a) is the best-fit model of equation (A9), while the dash-dotted 
and dashed lines in (b) are constant B and L lines derived by equations (All) and (A12) with the 
assumption of n c = 10 10 - 48 cm~ 3 . 




12 5 10 20 

<kT> (keV) 



Fig. 7. Plot of (kT) and (EM) in the flare phases. The symbols are the same as in figure 6. The 
dashed-dotted and dashed lines are constant B and L lines derived by equations (A15) and (A16) with 
the assumption n c = 10 10 48 cm~ 3 . 
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Table 1. Log of the Chandra ACIS-I Observations on the p Ophiuchi Cloud 



Obs. ID 


Sequence ID 


Date 


R.A.* 


Decl.* 


Exposure^ 








(J2000) 


(J2000) 


(ks) 


BF 


200060 


2000 Apr 13-14 


16 h 27 m 18?l 


-24°34'21"9 


100.6 


A 


200062 


2000 May 15-17 


16 h 26 m 35?3 


-24°23'12"9 


96.4 



The position of the detector aimpoint (the telescope optical axis). 
^The live time corrected from dcadtime. 
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Table 2. Chandra X-ray Sources in the p Oph Region. 



No. 


Count' R.A.T 


Decl.t 


(kTY 


\og{(EM)f 




Flux^ (i x ) s 


Red- 

x n 


Comment 




(J2000) 


(J2000) 


(kcV) 


(cm -3 ) 


(10 22 
cm -2 ) 







A-2* 



13.0(m)25 
v ; 55.22 



79359.2% 

03.02 



25 
37.3 



23 
36.1 



A-3 


829.9 


26 

07.05 


27 
24.4 


A-4 


23.7 


26 

07.36 


25 
31.4 


A K 
A-0 


OO/I 


9 (\ 

zo 

07.62 


9 7 

Z ( 

41.6 


A-0 


1 OQ 1 

iyo. i 


9 (\ 

ZO 

10.34 


on 

ZD 

54.6 


A-7 


83.1 


26 

12.42 


28 
49.0 


A-8 


16.2 


26 

13.30 


28 
23.0 


A-9 # 


12.9 


26 

14.34 


22 
28.1 


A-10« 


10.9 


26 

15.06 


25 
48.1 


A-ll 


59.8 


26 

15.71 


27 
49.3 


A-12 


180.3 


26 

15.81 


19 
22.2 


A-13 


16.2 


26 

16.32 


18 
43.0 


A-14 


2068.5 


26 

16.87 


22 
23.0 


A-15 


3741.7 


26 

17.06 


20 
21.6 



l.O(fixed) 


53.6 
56.3; 


51.9- 


45.5(>6.7) 


0.5 


38.4 


0.40(2) 




5.0(fixed) 


51.6 

53.2; 


50.7- 


15.0(1.6- 
129.3) 


0.6 


0.6 


0.31(2) 




2.9(2.8- 


54.8 


54.8- 


1.1(1.0-1.1) 


1257.2602.5 


1.51(436) 


Flare 1 


3.0) 


54.9; 














2.8(2.7- 


54.8' 


54.8- 




1016.0493.6 




Flare 2 


2.9) 


54.8; 














2.1(1.8- 


53.1 


53.0- 


3.7(3.2-4.2) 


13.7 


13.6 


0.81(76) 


Quiescent 


2.7) 


53.3; 














2.5(1.8- 


53.3 


53.1- 




26.1 


22.4 




Flare 


4.3) 


53.5; 














19.4(>0.2) 


51.4 
59.7; 


51.0- 


4.6(1.9- 
44.5) 


0.8 


0.4 


0.34(7) 




2.1(1.8- 


53.1 


53.0- 


2.8(2.5-3.1) 


13.4 


11.7 


0.72(114) 


Quiescent 


2.4) 


53.2; 














3.5(2.5- 


53.2 


53.1- 




34.5 


20.8 




Flare 1 


4.9) 


53.4; 














2.6(1.9- 


53.2 


53.0- 




22.0 


16.0 




Flare 2 


4.0) 


53.3; 














2.3(1.4- 


53.0 


52.6- 


5.8(4.2-8.5) 


8.2 


9.5 


0.74(17) 




4.3) 


53.4; 














3.3(1.4- 


52.1 


51.8- 


5.5(3.5- 


2.0 


1.6 


0.57(24) 




10.3) 


52.8; 




10.2) 










2.0(>0.2) 


51.4 

56.i; 


50.6- 


3.0(0.8- 
10.2) 


0.3 


0.3 


0.95(2) 




l.O(fixed) 


53.0 
54.1 


51.4- 


26.6(3.4- 
92.8) 


0.2 


8.9 


0.57(2) 




5.0(fixed) 


51.4 
51.9; 


50.4- 


10.8(0.8- 
46.5) 


0.4 


0.3 


0.39(2) 





3.7(>0.7) 


52.0(51.6- 


6.6(3.3- 


1.5 


1.3 




54.2) 


25.6) 






1.2(1.1- 


52.6(52.3- 


2.6(2.0-3.0) 


1.8 


3.1 


1.6) 


52.9) 








4.3(>0.9) 


52.8(52.4- 




15.1 


8.1 




53.3) 








5.2(>0.5) 


51.4(50.9- 


5.8(1.6- 


0.5 


0.3 




55.4) 


43.6) 






1.7(1.5- 


53.3(53.2- 


1.9(1.8-2.1) 


18.6 


16.3 


1.9) 


53.3) 








2.6(2.2- 


53.5(53.4- 




54.0 


33.4 


3.1) 


53.6) 








1.3(1.2- 


53.1(53.1- 


0.5(0.5-0.6) 


19.0 


10.7 


1.3) 


53.1) 








1.4(1.3- 


53.1(53.1- 




22.5 


12.1 


1.5) 


53.2) 








1.5(1.3- 


53.4(53.3- 




41.4 


20.7 


1.7) 


53.4) 









0.56(14) 

0.73(15) Quiescent 
Flare 

0.01(1) 

1.14(165) Quiescent 
Flare 

1.29(181) Quiescent 
Flare 1 
Flare 2 
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Table 2. (Continued) 



No. 



Count" R.A.T DecD {kT^ 



log((EM)y 



FluxS (i x ) s Red 
v 2 H 



Comment 



(J2000) (J2000) (kcV) 




(10 22 
cm -2 ) 



A-16 



126.2 26 



17.13 



12 
38.9 



1.5(1.4 
1.6) 
0.9(0.6 
1-2) 



53.2(53.2 
53.3) 
52.4(52.2- 
52.8) 



1.8(1.4-2.2) 1.1 2.5 



28.4 14.3 



1.26(16) 



Flare 3 



Background-subtracted X-ray counts in 0.5-2.0 keV, 2.0-9.0 keV, and 0.5— 9.0keV for soft-band, hard-band, and the other sources, 
respectively, (m) denotes sources with marginal detections (the confidence level < 99.9%, see subsection 3.1). Although the 
confidence levels of A-48 and A-H2 are significant enough, we regard them as marginal sources because of the larger source size 
(A-48) and severe contamination from A-2 (A-H2). 
' Right ascension and declination for all sources are 16 h and -24°. 
t Parentheses indicate the 90% confidence limits. 

^Observed flux (10 -14 erg s _1 cm -2 ) and absorption-corrected X-ray luminosity (10 29 erg s — 1 ) in 0.5-9.0 keV. 
II Rcduced-x 2 for the spectral fittings. Parentheses indicate the degrees of freedom. 

#Wc determine spectral parameters with fixed temperatures of 1 keV and 5 keV (see text). For sources which only show the 
parameters for kT = 1 or 5 keV, no good fitting is obtained for the other temperature. 

We assume the same abundances as the "F2" phase in Imanishi et al. (2002). The Quiescent spectrum is not obtained because 
the decay phases of the two flares occupy all of the light curve. 
t^The pileup effect is not corrected. 
fl-No spectral fit is done due to the limited statistics. 

^Abundances are free parameters (see subsubscction 4.8.1 in Paper I and Imanishi 2003). 

II II We make the flare spectra with a bit larger timescale in order to obtain as good statistics as possible. Errors of (EM) for A-29 
in the quiescent are not determined because of the limited statistics. 

**The spectra show possible edge absorption of neutral Ca or warm Ar. The nonthermal model also well reproduces the spectra 
(Hamaguchi, Imanishi 2002). 

***We assume the same temperature because of the limited statistics, 
t' 'These show non- thermal spectra (Imanishi 2003). 

* » * Imanishi ct al. (2002) proposed two temperature model with an unusual abundance pattern. 
§§§ Foreground star. The distance is 60 pc (Festin 1998). 

II II II The best-fit value of (kT) is not determined (larger than 10 keV), hence we assume 10 keV temperature for the estimation of 
the other parameters. 

###7Y H and reduced-x 2 are estimated by the simultaneous fittings with the identical sources in obs.-A. 
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Table 3. Identifications of the X-ray sources. 



No. 


Offset* 
(") 


Radio 1, 


— X-ray — -* 
PSPC HRI ASCATKT 


Other namcs§ 


Classll 


A-l 












A-2 


0.33 


R8, S7 


13 All 3 


DoAr21, ROXs8, YLW26, 

vjrOOZO, ±!jliaS14, orvol-D, loVJIU 


III 




n 4n 

U.^rU 






louio 


TT 
11 


A A 




xw ; 








A 


n fin 




U? A 1 Q 




TTT 
111 


A fi 

ri-D 


U.ZO 




ID 


pqcoc PRDtji: QI<"Q1 R TC!ni7 


TT 
11 


A 7 












A 8 












A-9 






... 




... 


A-10 












A-ll 












A-12 


0.47 






ORBR9 SKS1-7 IS018 


Ill 


A-13 












A- 14 


0.11 


Lpl 


18 A14 


ROXsll, GSS29, Eliasl8, 
CRBR10, SKS1-8, IS019 


II 


A-15 


0.44 




17 A15 


DoAr24. ROXslOA, GSS28, 
Eliasl9, CRBR11, SKS1-9, 
ISO20 


II 


A-16 


1.08 






2MASSI J1626172-241238 


? 



Offset between the Chandra and nearest 2MASS sources. For sources having a radio counterpart only, the offset between the 
Chandra and radio sources (Andre et al. 1987; Leous et al. 1991) are shown, which is indicated by (R). 

tRadio sources. S, R, and L indicate sources listed in Andre et al. (1987) (ROC), Stine et al. (1988) (SFAM), and Leous et al. 
(1991) (LFAM), respectively. Lp denotes possible radio sources in LFAM. 

*X-ray sources detected with ROSAT/PSPC (Casanova et al. 1995), ROSAT/HRI (Grosso et al. 2000), and ASCA (Kamata et al. 
1997). "IKT" denotes Chandra sources already listed in Paper I. 

^Source names used in the literature. Abbreviations for names are SR (Struve, Rudkjobing 1949), DoAr (Dolidze, Arakelyan 
1959), ROXs (Bouvier, Appenzeller 1992), S (Grasdalen et al. 1973), YLW (Young et al. 1986), GY (Greene, Young 1992), WL 
(Wilking, Lada 1983), VSSG (Vrba et al. 1975), GSS (Grasdalen et al. 1973), Elias (Elias 1978), BBRCG (Barsony et al. 1989), 
CRBR (Comeron et al. 1993), SKS (Strom et al. 1995), IRS (Wilking et al. 1989), and ISO (Bontemps et al. 2001). For sources 
having NIR counterpart only in the 2MASS catalog, we alternatively show the 2MASS source names. 

'Source classifications; I: class I, II: class II, III: class III, III C : class III candidate (table 5 in Bontemps et al. 2001), BD: brown 
dwarf, BD C : brown dwarf candidate (Imanishi et al. 2001b), F: foreground star (Festin 1998). Unclassified NIR sources are 
indicated by "?", and sources with no NIR counterpart are called "unidentified sources", labeled no data point (...) in this column. 
For BD/BDcS, we show the available IR classification in Bontemps et al. (2001) in the parentheses. 
#Offset from a nearest NIR source in Greene and Young (1992). 
"Grosso (2001). 
ttTsuboi et al. (2000). 

ttThe position of GY5 (Greene, Young 1992) is slightly shifted (~2."3) from the 2MASS source. 
§§A candidate of HH object (Gomez et al. 1998). 

I II From NIR spectroscopy, Wilking et al. (1999) derived for these M dwarfs masses higher than the hydrogen burning limit. 
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Table 4. Detected X-ray Flares. 



No 

1NU, 




7"r 


T d 


n 












I L.U LUIL J\.o J 










Place T 








"iX 

•Jo 


u.y 


199/09 icnl 
iz.z^y.z— 10.U ) 


u.o^u.o— l.U J" 




RF 9fi 
I J 1 -ZD 


1 X 


7 "iff; n o n"i 
< . o^o.u y.u ) 


0.0^4.4 O. / J 


7 n/'fi 4 7 ^\ 
1 .yj[\j.'± 1 .0) 




RF ^1 
or -01 


oZ 


9 0/9 c qol 
z.y (^z.o— o.z ) 


14.U^lz./— 1D.0J 


1 U . ( ^ 1 u . u~ 1 1 . ) 




RF £/l 
I ) 1 -04 


99 
ZZ 




1 /I «/1 °. 9 1 fi 1 "1 
14. 0(^ lo. Z— 10. 1 ) 


ii.4^y. < — lo.u j 


r idrc 1 


t 


07 0/07 1 07 

( .6(0 1 .1-0 I .DJ 3 




97 7/9/1 Q 01 c\ 

Z 1 .1 yZ'i.o— ol.D ) 




r lai c z 


RF SO 
or -oy 


1 

•i 


l.Zl U.O l.O J 


9 /I/O Q ,( 11 

o.4^z .y 4. 1 j 


O.Z^O.O O. 1 J 


Fl-iro 1 

r idi c 1 




GO 


4:.U^Z.4 O.Z J 


Q 0/9 1 7 c~\ 




FIqt-ci 9 
r idi c z 




t 


Z.o^l.O— 0.4 J 






r idic 








Place tt 

^/IclSS 11 








Z4: 


0. 1 ^4.G O.O J 


q 7/1 8 Q 
0. / ^1 .0 y.o ) 






A 1 A 


x 




A c/q q c 7"\ 
4.0^0.0 0. / J 


1 9 a/1 n fi 1 v^i 

IZ.o^lU.O 10. / J 


1 "i . ^ 1 . 1 . j 




A 1 ^ 
A- ID 


Z 


i 


1 1 7//I c 90 9~\ 
11. / ^4.0— Zo.Z J 


Zc5. 1 ^ZZ.l _ oU.o J 


r idic 1 




ou 


Q AO Q C 9"\ 


a 1 /q n ^ s"\ 




r idi c z 




8fi 


O.D^O. / — 1Z.1 ) 


1 c 1 /o 1 /1 r\\ 
lO.l^o.l— 4o.U J 




r idic 


A-zz 


9(1 


A 9/9 8 R 1 ^ 
4. Z^Z.Q— 0.1 J 


lU.o^O. O-10.0 ) 


/ .0(^0. /—o.z ) 




A 9^ 
A-Zo 


90 


n x* 

U.o 


4.U^O.O— 4.0 ) 


9 fl/1 7 9 ^ 
Z.U^l. / — Z.o J 


r id,i c 1 




7 1 

/ 4 


1 oYn s ill 
i.y^u.o 0.4 j 


Q 0/9 1 c o\ 
o.o^z.l 0.0 J 




r idi c z 


A 9/1 


79 
( Z 


9 QI'O 7 Ol^ 

z.y (^z. /— 0. 1 j 


190/11 q iqq"\ 
IZ.o^ll.o— lo.oj 


o.4^o.y— / .uj 




A 9^ 


1 (1 
±U 


n Q* 
u.y 


1 -O 


1 ( . u ^ 1 . y 1 . u j 


Fl-iro 1 
r idi c 1 




QQ 
•JO 


7.5(4.2-12.3) 


3.0(0.7-10.1) 




r lai C Z 




XX 

00 


1.1(0.5-1.7) 


5.6(3.8-9.1) 




FIqt-ci ^ 
r idi c 


A 

JA-<) I 


•JO 


3.0(1.5-5.8) 


8.1(5.6-11.8) 


z.^z.u-z.o J 




A 


89 

oz 


2.2(1.2-3.6) 


5.1(2.8-10.0) 


1 n fiM no 1 1 ^ 
1 u . ^ 1 u . u 1 1 . ) 




A fiQ 


7(1 
1 u 


26.9(11.0-121.3) 4.7* 


o.u 




A-/ / 


ox 

yo 


3.5(2.2-5.5) 


.J 


n 7/n p; n c^ 1 
U. / (^U.O— U.o J" 




A 78 

jtV- ( O 


■ J 4: 


1.5* 


3.1* 


1 Off! R-1 1 ^ 

l.U^U.O 1.1 J 




A 7Q 
a- ( y 


19 


2.8(2.5-3.0) 


9.0(20.8-16.3) 


9 zlYl 7 ^ 1 ^ 


r idi c 1 




R9 
OZ 


8.1(6.9-9.2) 


17.3(35.9-33.4) 




r idi c z 


RF S 
I > I -o 


00 


1.0* 


1.3* 


U.O^U.O U.O ) " 




BF-17 


32 


1.4(0.8-2.1) 


1.7(0.8-2.9) 


1.0(0.8-1.1) 




BF-27 


86 


1.5(0.6-2.5) 


2.1(1.1-3.4) 


0.9(0.7-1.1) 




BF-28 


18 


0.9(0.2-1.3) 


2.7(1.9-4.0) 


1.1(0.9-1.3) 




BF-31 


54 


2.5(2.3-2.7) 


9.1(8.4-10.0) 


1.1(0.8-1.3) 
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Table 4. (Continued) 



1N0. 


4- 

(ks) 


(ks) 


(ks) 


Q 

(count ks — 1 ) 


Comment 


Dr -OO 




hyp; o £ o^ 

0.U(0.Z— O.y J 


10.o(lo. ( — lo.z ) 


Z.o^Z.o— o.Z ) 




Dr -4Z 


i 


( .0(0.0— 10.1 ) 


A A 1 4 7 7 fl\ 
4.4(1 . 1—I.O) 


Q Q/Q A A 0\ 
0.0^0.4— 4. Z ) 




DT? CI 

Dr -01 





i n* 

l.U 


i n* 

l.U 


U.0^U.4— U.o ) 11 


r lare i 




p. 1 
01 


i ^* 

1.0 


1 -i* 

1 . o 




ri 01 .n 

r lare z 


or -oy 


oo 


Q O/'O 7 Q £^ 
O.Z(Z. |— O.O J 


1 .0(0. z— y.o J 


Z.\J[1. /— Z.o J 




or -Do 


04 


i 


o.y (o.o— o.y j 


Z.y ^Z.O— o.Z ) 




Dr -00 


yo 


o.U(l.z— 0.8 J 


u.o 


O /I/O 1 O 7\ 

Z.4(^Z.l— Z. ( J 




rSr - ( a 


4 


Q 7f0 K K fi\ 
o. 1 (z. 0—0.0 j 


7 Of£ Q O 

( .y (o.o— y.z ) 


zi.o^zu.y- zz.o ) 




BF-87 


84 


3.4(3.1-3.8) 


6.0(5.4-6.7) 


1.0(0.8-1.2) 




B.b-88 


14 


21.5(10.7-56.0) 


7.0(3.6-15.6) 


26.0(14.7-32.7) 


T?l „ 1 
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40 
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1 n/n 7 1 0\ 
1 .U(U. 1 l.z J 




IDT? (\Ci 

tit -yo 


•_> o 

oz 


59.2(41.9-96.7) 


13.8(10.8-17.9) 


07 o/O/i on o\ 
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— Class III C — 






A 1 O 

A-lz 


on 

20 


0.9* 


0.9* 


1 0/1 1 1 

1.3(1.1-1.5 j 




TDT? TO 

Dr - / Z 




..} 


4.1(2.1-6.8) 


1.2(1.0-1.4] 








- Unclassified NIR sources — 






A-26 


2 


3.1(1.0-57.9) 


4.0(3.1-5.1) 


1.4(1.0-1.7) 


Flare 1 




18 


1.3* 


4.6(2.3-10.9) 




Flare 2 




76 




20.1* 




Flare 3 


A-28 


50 


8.5(4.3-17.8) 


23.8(9.1-55.4) 


0.8(0.3-1.1) 
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Table 4. (Continued) 



1NO. 










Comment 




(tCSj 


(ksj 


(KB J 


(count ks ) 




A-7o 


TO 
( Z 


1.9* 


24.0(12.2-129.1) 


1.0(0.o-1.2j 




nil r r 


•3 (J 


1.5* 


0.8* 


0.7(0.5-0.9)" 




r>r -Oz 


on 


3.7(2.6-5.0) 


6.9(4.9-9.9) 


1 .z(U.o— l.o J 


r lare i 




Do 


5.7(5.0-6.5) 


10.5(9.3-12.0) 




r idre z 






- Unidentified sources — 






A-29 


80 


1.0* 


2.1* 


0.5*11 




BF-36 


30 


4.4(1.7-9.4) 


.J 


0.8(0.6-0.9) 




BF-92 


20 


0.5* 


2.6(1.1-4.4) 


1.9* 





We can not determine errors because of the limited statistics and/or strong coupling of the other parameters. 
^Fitting is obtained for light curves in 2."5— 12."5 (A-2) and 2."5-7."5 (BF-64) radius circles in order to avoid 
the pileup effect. 

^No fitting is possible because of the limited statistics and/or the boundary of the observations. 

§We relax the flare peak time (t p ) to be free in order to obtain as good results as possible. 

II This should be the upper limit because the fitting does not include the time bin with zero counts. 



Table 5. Results of two-sample tests for the differences of the observed parameters among classes. 



Parameter — Sample size*— - Probability for the null hypothesis^ — 

I II III+III C I VS II+III+IIIc I+II VS III+III C II VS III+III C 











GW 


logrank GW 


logrank GW logrank 


logix (Flare) 


8 


36 


16 


0.059 


0.063 




logix (Quiescent) 


8 


61 


26 








(kT) (Flare) 


7 


35 


16 


0.079 


0.060 


0.017 ... 0.047 


(kT) (Quiescent) 


7 


61 


24 


0.008 


0.022 0.046 


0.037 ... 0.026 


T~r 


6 


28 


8 










7 


28 


12 









The number of data points used for the two-sample tests. 
T Probability that the hypothesis of two distributions being the same is true, which is derived with the Gehan's 
generalized Wilcoxon test (GW) and the logrank test (Feigelson, Nelson 1985). Blank (...) indicates that the 
probability is larger than 0.1. 
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Table 6. Estimated Mean Values of the Flare Physical Parameters. 



Class 




This work 




— Kl VS 


hjM - 




iog(n c ) l + 


iog{By 


log(L) 1 


iog(issY) 3 


lOg(LsY) 3 




(cm -3 ) 


(G) 


(cm) 


(G) 


(cm) 


I 


10.48±0.09 


2.7±0.1 


10.3±0.1 


2.6±0.1 


10.4±0.2 


TT 

11 




Z.o±U.l 


1U.4±U.1 


Z.OztU.l 


1 C\ AJ-f\ 1 

lU.4±U.l 


III+III C 




2.3±0.1 


10.4±0.2 


2.2±0.1 


10.7±0.1 






(2.4±0.1)H 


(10.0±0.1)H 






Equation* 


(5) 


(All) 


(A12) 


(A15) 


(A16) 



Equations in the text used for the estimation of each parameter. 
^The values are derived by assuming Ma = 0-01 
+We assume the same values of n c for all classes (see text). 

§These values are estimated using the derived values of n c (= 10 10 - 48 cm -3 , column 2). 

II The mean values for class III+III C when we exclude the flares with unusually long 

timescales (A-2, A-63, and BF-96). 
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